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Isocurvature perturbations naturally occur in models ol inflation consisting of more than one scalar 
field. In this paper we calculate the spectrum of isocurvature perturbations generated at the end of 
inflation for three different inflationary models consisting of two canonical scalar flelds. The amount 
of non-adiabatic pressure present at the end of inflation can have observational consequences through 
the generation of vorticity and subsequently the sourcing of B-mode polarisation. We compare two 
different deflnitions of isocurvature perturbations and show how these quantities evolve in different 
ways during inflation. Our results are calculated using the open source Pyflation numerical package 
which is available to download. 
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I. INTRODUCTION 

Inflationary models consisting of more than one scalar 
field have been long known to create isocurvature, or 
non-adiabatic pressure, perturbations due to the relative 
entropy perturbation betvifeen the different fields [J [5] . 
Isocurvature perturbations during inflation are interest- 
ing since, for example, they induce the evolution of the 
curvature perturbation on super-horizon scales [3H5] , and 
can give rise to large non-gaussianities [51 E] ■ 

The fraction of isocurvature perturbations present 
around recombination is derived from observations of the 
cosmic microvifave background (CMB) such as those made 
by the Wmap satellite, and there have been several recent 
works studying this [SHTT]. However, even in the 'era of 
precision cosmology' the Wmap seven year results allow 
for some isocurvature, constraining the power spectrum 
of entropy perturbations to be of order 10% that of the 
usual adiabatic perturbations |12| . 

The most common method for computing perturba- 
tions from inflation is to use the so-called '(57V formal- 
ism' which is a gradient expansion relating the 
perturbed number of e-folds to the curvature perturba- 
tion on uniform density hypersurfaces. This is a concise 
method, but has some draw-backs, primarily that it is 
only valid on super-horizon scales. 

An alternative method for computing inflationary per- 
turbations is to use the full cosmological perturbation 
theory, a powerful technique which has been honed dur- 
ing the past few decades [TMIS] (see also Refs. [El ED] 
and references therein for a comprehensive list). It is this 
approach which we adopt in this article. This will enable 
us to consider the non-adiabatic pressure perturbation 
during general multi-field inflation and obtain its power 
spectrum at the end of inflation. 



There are several reasons for wanting to consider the 
non-adiabatic pressure perturbation produced towards 
the end of inflation in such a comprehensive manner. 
One interesting reason is the recent realisation that a 
non-adiabatic pressure perturbation can source vortic- 
ity at second order in cosmological perturbation theory 
pm23j . This vorticity will, in turn, have an effect on 
the CMB, since vector perturbations naturally source B- 
mode, or rotational, polarisation of the CMB radiation 
|24) . They can also possibly source primordial magnetic 
fields '251 thus providing another way in which to dis- 
tinguish between the many different inflationary models 
and help to constrain the inflationary paradigm. 

In this article we will consider several models of infla- 
tion consisting of two canonical scalar flelds which pro- 
duce a non-negligible non-adiabatic signal, and calculate 
the spectrum of entropy perturbations for these mod- 
els in two different ways. The first quantity is the non- 
adiabatic pressure perturbation defined in terms of the 
total fluid perturbations of pressure and energy density 
and is denoted ^Pnad- This is compared with the spec- 
trum of the rotated fleld values in the isocurvature direc- 
tion Ss. We consider our results in comparison with those 
previously obtained and comment on the how well the 
models we investigate are constrained by current obser- 
vational data. We also show that the full non-adiabatic 
pressure perturbation, dPnad, evolves in a quite different 
way to Ss. It is important to distinguish between these 
quantities in evaluating the potential of particular phys- 
ical models when considering the generation of vorticity. 
The numerical code used in these calculations has been 
released under an open source license and is available to 
download [26 . 

This paper is organised as follows: in the next section 
we derive expressions for the non-adiabatic pressure per- 
turbations in multi-fleld inflation. We then outline our 
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numerical method in Section III and present our results 
in Section llVl We discuss our conclusions and directions 
of possible future work in Section jV} 
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II. PRELIMINARIES 



and total adiabatic sound speed as 



An adiabatic system is one in which 
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for any scalars X and Y , and where an overdot denotes 
the derivative with respect to coordinate time, t. It is 
usual to describe a fluid using its pressure and energy 
density, P and p, and so an adiabatic fluid is defined as 
one in which 



Thus, we arrive at the following expression for the non- 
adiabatic pressure perturbation in the system consisting 
of the two fields Lp and x 
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In general, though, a fluid system is not adiabatic, and 
so the pressure perturbation can be expanded as 



Using one of the Einstein fleld equations we can relate 
the lapse function (j) to field variables as 



,5P = (5P„ad+C^,5p, 



(2.3) 



H(j) = ■iTTG{ip5ip + x5x) ■ 



where JPnad is the non-adiabatic pressure perturbation a-nd so Eq. ( 2.11 ) becomes 
and Cg = P I p is the adiabatic sound speed for the fiuid. 



For a system consisting of more than one fluid or field, 
this can then be expanded further into an intrinsic and 
relative entropy perturbation |17j 
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In this article we do not consider this split, and instead 



we use Eq. (2.3) to calculate the value of <5Pnad as 



nad 



5P~cl5p. 



(2.5) 



It is this definition of the non-adiabatic pressure per- 
turbation that we calculate in the next section. In order 
to compare with the comoving curvature perturbation TZ, 
which is (in the flat gauge) 



In order to calculate the pressure and density of the 
system we consider linear, scalar perturbations to a flat 
Friedmann-Robertson- Walker spacetime, and choose the 
gauge in which spatial hypersurfaces remain unperturbed 
- the uniform curvature gauge. The line element thus 
takes the form 

ds^ ^ -(1 + 2(l))df + 2a{t)B ,dtdx' + a'^{t)5^jdx'dx^ , 

(2.6) 

where a{t) is the scale factor, <j> is the lapse function 
and B is the scalar shear. We consider a system of two 
canonical scalar fields with the Lagrangian 



n 



(2.14) 



we use a comoving entropy perturbation S introduced in 
Refs. EHl and defined as 



S = -pSPnad 



(2.15) 



C 



1{^' + X')+V{^,X)- (2.7) 



An alternative way to calculate perturbations in two- 
field infiationary models is to perform the field rotation 
into an adiabatic field, a and an isocurvature field s as 
in Ref. 128J: 



The total energy density and pressure perturbations 
can be written in terms of the fields as [191 l27l 



and 



6(T — cos 9S(p + sin 66x , 
6s = — sin 65ip + cos 0Sx , 



(2.16) 
(2.17) 



Sp^^(ipaSip^~ ipl(t>+VaSipa) , (2 



where tan(0) — x/ip and a comoving isocurvature per- 
turbation can be then be defined as 



S=-Ss, 
a 



(2.18) 



<5P = ^ (^PJ^^ - y'o^ - y.gS^^ , (2.9) 



where & = \/ (p^ + X^i The two isocurvature perturba- 
tion definitions S and S are then related to one another 



3 



through dH] ^ 
S = 



■ 8ttG . 2 r 

d(T a da 

2H 



aba 



(2,Ea^e + 2aV^„B)bs \ , 



(2.19) 



where V^a — 'Ylia'^-OL^al ^ and we have substituted for 
= ^-nGaba. In Section IV results are given for both 
5 and 5 to show the differences in evolution. In the 
slow roll and large scale limits S and S become equiva- 
lent and bs is an isocurvature direction. In general this 
is not the case as can be seen in Section IIVI when slow 
roll breaks down at the end of inflation. It is clear from 



Eq. (2.19) that even when the isocurvature perturbation 



bs is highly suppressed, S can still evolve due to the de- 
pendence on the adiabatic component ba. This is an im- 
portant distinction between the two definitions of isocur- 
vature which is not often recognized. 



III. NUMERICAL PROCEDURE 

We have used the Pyflation numerical package to 
obtain the results in this paper (2^1 HO]. This package 
uses the Numerical and Scientific Python libraries [3TH33] 
and is based around a Runge-Kutta ODE solver [M] . The 
Pyflation package has been updated to solve first order 
systems of multiple inflationary fields and this new ver- 
sion is now available for download under an open source 
license [55]. 

As described above we evolve the perturbed field val- 
ues bif> and bx in contrast with other approaches. In fact 
the quantum nature of the perturbations on sub-horizon 
scales means that for a multi-field system the operators 
for each of the N fields contain N independent annihila- 
tion operators as 



(3.1) 



Our numerical system solves the evolution equation for 
the N X N mode function matrix ^a/s following Ref. [3S] . 
It is common to replicate this behaviour without em- 
ploying the full mode matrix by running the simulation 
multiple times, with each of the fields in turn set to zero 
initially. If the matrix approach or the multiple runs ap- 
proach are not used the numerical results will be incom- 
plete as they will not account correctly for the cross-terms 
in the mode matrix (e.g. ^0,^, in the two-field case) (see 
Refs. [36l |37j and the notes contained there on previous 
numerical calculations). 



^ Note that the final term in this expression differs sUghtly from 
that in Eq. (42) of Ref. |28| . In the slow roll limit when a is 
negligible the expressions are equal. 



The Klein-Gordon equations for the scalar field pertur- 
bations must be rewritten in terms of the mode function 
matrices [3S]. For clarity these equations are presented 
here using coordinate time, however in the numerical sys- 
tem the time variable is the number of efolds A/", and all 
the ODEs are expressed in terms of derivatives with re- 
spect to Af. For the first order mode function matrix 
elements the Klein-Gordon equation is 



^al3 + iH^af} 



SttG { . . 8nG . 



= 0. 



(3.2) 



The initial quantum state is taken to be the Bunch- 
Davies vacuum state with the scalar fields uncorrelated 



aV^k 



(3.3) 



where 77 = J dt/ a is conformal time. All the quantities 
calculated from the field perturbations inherit this mode 
matrix structure so for example bPnud = J2a ^Pnada&a- 
When using the mode matrix approach it is important 
to calculate the power spectra of quantities correctly us- 
ing the commutation relations of the quantum operators 
in order to account for the cross-terms. For example the 
power spectrum of comoving curvature perturbations V-jz 
is given by 



2n' 



(3.4) 



IV. RESULTS 

In this section we present our numerical results on 
three inflationary potentials involving two scalar fields. 
These potentials have been well studied in the literature 
and we have used these previous results to ensure our 
code is working correctly. 



A. Double quadratic inflation 



We first consider 
quadratic inflation [51 
form 



the well-studied case of double 
] for which the potential takes the 
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(4.1) 



We consider the particular case when m-^ — 7m^ as stud- 



ied in Refs. [Ml [37], and we set = 1.395 x IO^^A/pl 
in order for the curvature power spectrum at the pivot 
scale to match the WMAP7 maximum likelihood value 
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FIG. 1: A comparison of the power spectra Vn (red straight FIG. 3; A comparison of the power spectra of 6P (red straight 
line), Vs (green dashed line) and (blue dotted line) for line) and iJPnad (green dashed line) for the double quadratic 
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FIG. 2: A comparison of the power spectra Vn (red straight 
line), Vs (green dashed line) and (blue dotted line) for 
the double quadratic potential in terms of k at the end of 
inflation. 



of 2.45 X IQ-^. For the results shown here in Figures Tj2 
and [3] the initial field values are <^o = Xo = 12MpL and 
the derivatives of the fields are given their slow roll val- 
ues. The spectral index of the curvature perturbations 
for these parameter choices is approximately n-fi ~ 0.937 
(with no running allowed). 

In this model the inflationary dynamics are originally 
dominated by the x field until around 30 efolds before the 
end of inflation when the fleld becomes dominant. This 
is associated with a rise in the curvature perturbation 
amplitude as shown in Figure [l] where the power spectra 
for the comoving curvature perturbation and the two def- 
initions of isocurvature perturbations are compared. In 



all the graphs in this section results are plotted against 
either the number of efoldings A/" left until the end of 
inflation (Af) or the comoving wavenumber k of the per- 
turbation modes in units of inverse Megaparsecs. When 
plotted against J\f the mode shown is the WMAP pivot 
scale k = 0.002Mpc-^ 

The evolution of S matches that found in Ref. [37] 
with the amplitude dropping off signiflcantly after the 
change-over. This is a good consistency check for our 
mode function matrix approach as that paper used the al- 
ternative multiple-run method as described in Section III 
and performed the field redefinition into a and s before 
the numerical run. ^ 

The evolution of S is not similar to S especially as the 
mode crosses the horizon around 60 efoldings before the 
end of inflation. The amplitude of S continues to reduce 
for a few efoldings after horizon crossing before rapidly 
increasing to reach a peak at the cross-over time. It then 
rapidly drops again and at the end of inflation is many 
orders of magnitude smaller than TZ as shown in Figure [2] 
The amplitude of S is even further suppressed in com- 
parison but the numerical value at such a relatively small 
level must be treated with caution. We have also calcu- 
lated the spectral index of the isocurvature perturbations 
using a similar expression to n-ji: 

alog(fc) 

As one might guess from Figure[2]the two spectral indices 
have very similar values. This might be expected given 
the dependence of SP^ad and therefore S on the com- 
binations of the fields which are also present in TZ, but 
this result will be very useful for future studies of the 
generation of vorticity as discussed further in Sectio n W] 
The quantity of interest we described in Section 11] is 
the non-adiabatic pressure perturbation 5Pnad and this 
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FIG. 4: A comparison of the power spectra Vu (red straight 
hne), Vs (green dashed line) and Vq (blue dotted line) for 
the double quartic potential at the Wmap pivot scale. 



FIG. 6: A comparison of the power spectra V-r. (red straight 
line), Vs (green dashed line) and V§ (blue dotted line) for 
the double quartic potential at the Wmap pivot scale over 
the last 5 e-foldings of inflation. 
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FIG. 5: A comparison of the power spectra Vn (red straight 
line), Vs (green dashed line) and V^ (blue dotted line) for the 
double quartic potential in terms of k at the end of inflation. 
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FIG. 7; A comparison of the power spectra of SP (red straight 
hne) and SPnud (green dashed line) for the double quartic 
potential at the Wmap pivot scale. 



is shown in Figure [3] along with the total pressure per- 
turbation SP. The evolution of (5Pnad is broadly similar 
to S, after a sharp fall on crossing the horizon the am- 
plitude increases until the cross-over point at which time 
it sharply falls again. This is to be expected as the rate 
of change in TZ is proportional to SPnad [19l I29j. 



B. Double quartic inflation 

Next, we study a special case of hybrid inflation that 
was considered in Refs. [371 133] • The potential for this 



model takes the form 



F(^,X)=A4 



^2 



(4.3) 



with the parameter values v — O.lOAfpL, (fie — O.OlA/pL 
and ^ = IG'^AfpL- In this case A can be normalised 
to match the WMAP results by setting A = 2.36 x 
10~'^Mpi^. The fields are started at the initial values 
ifo = O.OlMpL and xo = 1-63 x IQ-^AfpL. For these 
choices the spectral index at the end of inflation is around 
n-jz = 0.932 when no running is allowed. 

As shown in Figure |4] the behaviour in this case is 
markedly different. TZ continues to evolve outside the 
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FIG. 8: A comparison of the power spectra Vu (red straight 
hne), Vs (green dashed line) and Vq (blue dotted line) for 
the product exponential potential at the Wmap pivot scale. 



FIG. 9: A comparison of the power spectra V-r. (red straight 
line), Vs (green dashed line) and (blue dotted line) for 
the product exponential potential in terms of k at the end of 
inflation. 



horizon and S has a greater amplitude than TZ for a few 
efoldings after horizon crossing, again matching the re- 
sults of Ref. [37] . The evolution of S outside the horizon 
is unremarkable until around the final five efoldings of in- 
flation when a sharp increase in amplitude occurs. Close 
inspection of these final efoldings in Figure [6] shows that 
the amplitude of S rebounds to reach significant levels at 
the end of infiation. From Figure [5] we can see that Vs 
is of the order of a few percent of Vn . 

The evolution of the pressure perturbations is also very 
different from the standard double quadratic case. Fig- 
ure [7] shows 5P increasing throughout the super-horizon 
evolution while the level of (5Pnad is relatively constant 
until the final few efoldings. 



C. Product exponential 
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Finally we investigate a quite particular potential 
which has appeared recently in the literature in con- 
nection with the production of non-Gaussianity [40l [41] . 
This is a two field model with the product separable form 

V = VaLp^e-^^' , (4.4) 

where we have set A = O.G5/AfpL and normalized the 
resuhs to the WMAP value by setting Vb = 5.37 x 
10~^'^AfpL. For the results shown here the initial field 
values are ipa — 18Mpl and xo — 0.001A/pl. 

This model has been used as an example of a prod- 
uct separable potential in which it is possible to generate 
a large negative amount of non-Gaussianity^. Ref. [ID] 



^ The sign convention used is that of the WMAP team. 



FIG. 10: A comparison of the power spectra of 5P (red 
straight line) and 5Pnad (green dashed line) for the product 
exponential potential at the Wmap pivot scale. 



followed the evolution using the 5N formalism until the 
slow roll approximation becomes invalid before the end 
of inflation and showed that /nl values of the order of 
—35 can be obtained during infiation. In Ref. f5T] the 
authors showed that extremely large oscillations of the 
value of /nl occur after slow roll breaks down and on 
into the reheating phase. Because of this and the pres- 
ence of isocurvature it is not clear whether the values of 
/nl during inflation can be related to the observed values. 
In Figure [8] the late time evolution of TZ can be clearly 
seen as inflation ends. What is also very clear from our 
results is that the isocurvature perturbations measured 
by S have significantly larger amplitudes than the adia- 
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batic perturbations. The magnitude of S is smaller than 
TZ but it is a sizeable fraction at the end of inflation as 
shown in Figure [9] and puts this model in conflict with 
the restrictions on isocurvature fraction from the CMB. 

It is also clear from Figure |9] that the scale depen- 
dence of the curvature perturbations is not close to the 
maximum likelihood value from the WMAP 7 year data 
release. At the end of inflation the spectral index of the 
curvature perturbations is around n-jz ~ 0.794 when no 
running is allowed'^. However, the spectral index changes 
considerably in the last few efoldings of inflation. The ef- 
fects of reheating can not be neglected for this model and 
it is possible that after reheating the value of n-ji will be 
at least closer to the current observational limits. This 
is highlighted in more detail in Ref. j42] . 

As is now familiar the amplitude of i5Pnad drops as the 
mode crosses the horizon and TZ stops decreasing. The 
amplitude of (5Pnad then rises through the super-horizon 
evolution before crossing through zero as the amplitude 
of TZ reaches its maximum, and finally reaching levels 
comparable with SP. 

In this section we have shown the results of our numeri- 
cal simulations for three different two field potentials. We 
have seen that the evolution of S and SPnad is markedly 
different to that of S and 6s. We have focussed for the 
most part on models in which a significant fraction of 
isocurvature is present near the end of inflation and these 
models will provide a good starting point for the gener- 
ation of vorticity in the post-inflationary era. We have 
also seen the third model, the product exponential po- 
tential, is in some tension with the observational limits 
on isocurvature fraction and also the spectral index. The 
effects of reheating in this model are beyond the scope of 
this work but there are suggestions, especially from the 
results of Ref. ^] that the changes during this phase 
would be worthy of further consideration. 



V. DISCUSSION 

In this paper we have studied the spectra of entropy, or 
isocurvature, perturbations in various inflationary mod- 
els, focussing on models consisting of two scalar fields. 
We have undertaken a fully numerical study without re- 
sorting to an expansion in slow roll parameters. We 
have shown that the non-adiabatic pressure perturba- 
tion SPnad and its related quantity S can evolve in a 
quite different way to the isocurvature mode Ss which is 
frequently used. 

The study of isocurvature perturbations has been very 
popular recently. Other work includes looking at isocur- 
vature perturbations from multi-field infiation with non- 
canonical kinetic terms [SB] HI], or using the field 



is again very similar to n-fi in this case. 



redefinition, first considered in Ref. [55], into adiabatic 
and entropic modes [37l |45] . 

We have studied three inflationary potentials in this 
paper and described the evolution of the isocurvature 
perturbations until the end of infiation. After inflation, of 
course, the universe will need to go through a phase of re- 
heating during which the scalar fields driving inflation are 
converted into the standard model particles. This phase 
will have an effect on the spectrum of the non-adiabatic 
perturbations that we have obtained during infiation, and 
could cause the resonant growth of the entropy pertur- 
bations (see, e.g., Ref. [IB]). However, since there is no 
agreed upon mechanism of reheating, we do not study 
that here, instead focussing on determining which infla- 
tionary models in their own right can generate sizeable 
entropy perturbations. 

As mentioned in the Introduction, one of the motiva- 
tions for undertaking this study is the recent work high- 
lighting the importance of entropy perturbations on the 
generation of vorticity at second order in cosmological 
perturbation theory. As shown in detail in Ref. [21] , sec- 
ond order vorticity, uj2ij evolves (during radiation domi- 
nation) according to 

UJ21J - HuJ2ij OC (5p,[-,(5Pnadl,i] , (5-1) 

where the subscripts denote the order in perturba- 
tion theory and the square brackets denote anti- 
symmetrisation. So, vorticity is sourced by the coupling 
between linear energy density and entropy gradients. The 
perturbation to the energy density is well-known, and 
forms the basis of much cosmological study. However, 
entropy perturbations are rarely studied in cosmology in 
this context and so the results presented here will allow 
us to advance this field of study. 

At least for the models studied in this paper, the scale 
dependence of <5Pnad is very similar to that of the adia- 
batic curvature perturbations. This could have been an- 
ticipated but up to now this relationship has not been 
assumed in the literature on vorticity generation and 
instead various ansatze have been made for how (5Pnad 
varies with k. We hope that the results in this paper 
will be able to inform future work on the generation of 
vorticity and other non-linear effects using non-adiabatic 
pressure perturbations. 

As we have seen, for some models it is expected that 
reheating could strongly affect both the adiabatic and 
isocurvature results. The next step is to implement a 
straightforward reheating method, with the inflaton de- 
caying into a small number of species. It has also been 
shown that non-adiabatic pressure can be generated dur- 
ing a radiation phase even when the initial conditions are 
purely adiabatic J47]. By adding in the results presented 
here for isocurvature from the end of inflation it will be 
possible to refine these simulations and allow better com- 
parisons with the observational data at later times in the 
universe's evolution. 

In conclusion we have shown the evolution of two dif- 
ferent types of isocurvature perturbation for a range of 
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two field models of inflation. Our numerical results in- 
dicate that the evolution of (SPnad can be quite different 
from the isocurvature mode Ss which is usually consid- 
ered in the literature. Our results are significant when 
considered in conjunction with the generation of vorticity 
by (5P„ad, and the effect this may have on the interpreta- 
tion of observational results from the CMB, in particular 
B-mode polarisation. 



Acknowledgements 

The authors are grateful to Karim Malik, David Mul- 
rync, Paul Saffin, David Seery and Lixin Xu for useful 
discussions and for comments on a previous version of 
this manuscript. IH is supported by the STFC under 
Grant ST/G002150/1, and AJC is funded by the Sir 
Norman Lockyer Fellowship of the Royal Astronomical 
Society. 



L. A. Kofman, Phys. Lett. B173, 400 (1986). 
A. D. Linde, Phys. Lett. B158, 375 (1985). [27 
J. Garcia-Bellido and D. Wands, Phys.Rev. D53, 5437 
(1996), astro-ph/9511029. [28 

G. Rigopoulos and E. Shellard, Phys.Rev. D68, 123518 
(2003), astro-ph/0306620. [29 
D. Wands, K. A. Malik, D. H. Lyth, and A. R. Liddle, 
Phys. Rev. D62, 043527 (2000), astro-ph/0003278. [30 
D. Langiois, F. Vernizzi, and D. Wands, JCAP 0812, 004 
(2008), 0809.4646. [31 
M. Kawasaki, K. Nakayama, T. Sekiguchi, T. Suyama, 

and F. Takahashi, JCAP 0811, 019 (2008), 0808.0009. [32 
R. Bean, J. Dunkley, and E. Pierpaoli, Phys.Rev. D74, 
063503 (2006), astro-ph/0606685. ' [H 

J. Garcia-Bellido, p. AHEP2003/041 (2004), astro- [33 
ph/0406488. 

M. Bucher, J. Dunkley, P. Ferreira, K. Moodley, and | 

C. Skordis, Phys.Rev.Lett. 93, 081301 (2004), astro- [34 
ph/0401417. 

D. Parkinson, S. Tsujikawa, B. A. Bassett, and L. Amen- 
dola, Phys.Rev. D71, 063524 (2005), astro- ph/0409071. 

E. Komatsu et al. (2010), 1001.4538. [35 
M. Sasaki and E. D. Stewart, Prog.Theor.Phys. 95, 71 
(1996), astro-ph/9507001. [36 
D. H. Lyth, K. A. Malik, and M. Sasaki, JCAP 0505, 

004 (2005), astro-ph/0411220. [37 

D. H. Lyth and D. Wands, Phys. Lett. B524, 5 (2002), 
hep-ph/0110002. [38 
J. M. Bardeen, Phys. Rev. D22, 1882 (1980). 

H. Kodama and M. Sasaki, Prog. Theor. Phys. Suppl. [39 
78, 1 (1984). 

V. F. Mukhanov, H. A. Feldman, and R. H. Branden- [40 
berger, Phys. Rept. 215, 203 (1992). 

K. A. Malik and D. Wands, Phys. Rept. 475, 1 (2009), [41 
0809.4944. 

K. A. Malik and D. R. Matravers, Class. Quant. Grav. [42 
25, 193001 (2008), 0804.3276. [43 
A. J. Christopherson, K. A. Malik, and D. R. Matravers, [44 
Phys. Rev. D79, 123523 (2009), 0904.0940. 

A. J. Christopherson and K. A. Malik, Class. Quant. Grav. [45 
28, 114004 (2011), 1010.4885. 

A. J. Christopherson, K. A. Malik, and D. R. Matravers, [46 
Phys.Rev. D83, 123512 (2011), 1008.4866. 
W. Hu and M. J. White, New Astron. 2, 323 (1997), 
astro-ph/9706147. [47 

E. Fenu, C. Pitrou, and R. Maartens, 
Mon.Not.RoyAstron.Soc. 414, 2354 (2011), 1012.2958. 
L Huston, Pyflation: Cosmological perturbations for 



Python, http : //pyflation. ianhuston. net' (2011-). 
J.-c. Hwang and H. Noh, Class. Quant. Grav. 19, 527 
(2002), astro-ph/0103244. 

C. Gordon, D. Wands, B. A. Bassett, and R. Maartens, 
Phys. Rev. D63, 023506 (2001), astro-ph/0009131. 

K. A. Malik and D. Wands, JCAP 0502, 007 (2005), 
astro-ph/0411703. 

L Huston and K. A. Malik, JCAP 1110, 029 (2011), 
1103.0912. 

T. Oliphant et al., NumPy: Scientific computing tools for 
Python, http://numpy.scipy.org/ (2001-). 
T. E. Oliphant, Computing in Science & Engineering 
9, 10 (2007), URL |http : //link. aip . org/link/?CSX/9/| 

wn 

E. Jones, T. Oliphant, P. Peterson, et al., SciPy: Open 
source scientific tools for Python, http: //www. scipy . 
_org/ (200 1-). 

M. Abramowitz and L A. Stegun, Handbook of Mathe- 
matical Functions with Formulas, Graphs, and Mathe- 
matical Tables (Dover, New York, 1964), ISBN 0-486- 
61272-4. 

D. S. Salopek, J. R. Bond, and J. M. Bardeen, Phys. Rev. 
D40, 1753 (1989). 

Z. Lalak, D. Langiois, S. Pokorski, and K. Turzynski, 
JCAP 0707, 014 (2007), 0704.0212. 

A. Avgoustidis, S. Cremonini, A.-C. Davis, R. H. Ribeiro, 

K. Turzynski, et al. (2011), 1110.4081. 

D. Langiois, Phys. Rev. D59, 123512 (1999), astro- 

ph/9906080. 

H. Kodama, K. Kohri, and K. Nakayama, 
Prog.Theor.Phys. 126, 331 (2011), 1102.5612. 
C. T. Byrnes, K.-Y. Choi, and L. M. Hall, JCAP 0810, 
008 (2008), 0807.1101. 

J. Elliston, D. J. Mulryne, D. Seery, and R. Tavakol, 

JCAP 1111, 005 (2011), 1106.2153. 

M. Dias and D. Seery (2011), 1111.6544. 

X. Gao, JCAP 1002, 019 (2010), 0908.4035. 

C. M. Peterson and M. Tegmark, Phys.Rev. D83, 023522 

(2011), 1005.4056. 

S. Tsujikawa, D. Parkinson, and B. A. Bassett, Phys.Rev. 

D67, 083516 (2003), astro-ph/02 10322. 

R. AUahverdi, R. Brandenberger, F.-Y. Cyr-Racine, and 

A. Mazumdar, Ann.Rev.Nucl.Part.Sci. 60, 27 (2010), 

1001.2600. 

L A. Brown, A. J. Christopherson, and K. A. Malik 
(2011), 1108.0639. 



